Medium modification of jet fragmentation in Au+Au collisions at ^/s^ = 200 GeV 
measured in direct photon-hadron correlations 
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The jet fragmentation function is measured witli direct plioton-hadron correlations in p+p and 
Au+Au collisions at ^/s]^ = 200 GeV. The pr of the photon is an excellent approximation to the 
initial pr of the jet and the ratio zt ~ Pt /Pt is used as a proxy for the jet fragmentation function. 
A statistical subtraction is used to extract the direct photon-hadron yields in Au+Au collisions 
while a photon isolation cut is applied in p+p. Iaa, the ratio of jet fragment yield in Au+Au to 
that in p+p, indicates modification of the jet fragmentation function. Suppression, most likely due 
to energy loss in the medium, is seen at high zt- The fragment yield at low zt is enhanced at large 
angles. Such a trend is expected from redistribution of the lost energy into increased production of 
low-momentum particles. 

PACS numbers: 25.75.Dw 



Experiments at the Relativistic Heavy Ion Collider 
have observed the formation of a Quark-Gluon Plasma, 
reported as a fundamentally new state of matter 
High momentum quarks and gluons (partons) lose en- 
ergy as they traverse this matter, resulting in the ob- 



served suppression of high transverse momentum (high 
Pt) hadrons in central heavy- ion collisions [^-Q. 

Direct photons, however, escape the medium unmod- 
ified [§], since they do not interact via the strong force. 
This makes them an ideal probe with which to cali- 
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brate the energy of an initial hard scattering. At lead- 
ing order, direct photons are produced via the quantum- 
chromodynamics analog of Compton scattering, g + g — > 
(7 + 7. In the limit of negligible initial transverse momen- 
tum, the final state quark and photon are emitted back- 
to-back in azimuth with the photon balancing the trans- 
verse momentum of the jet arising from the quark. Deter- 
mining the initial momentum of the parton is key to mea- 
suring the fragmentation function of the quark jet. This 
initial momentum is provided by the measured energy of 
the unmodified direct photon, via direct photon-hadron 
correlations A study of direct photon-hadron cor- 

relations in p+p collisions at .^s„„ = 200 GeV has pro- 
vided a measurement of the fragmentation function in 
agreement with measurements from e 
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'e~ collisions 

In Au-|-Au collisions contributions from next-to-leading- 
order processes and medium induced photon production 
are expected to be smah («10%) at high px [12]. 

Parton energy loss in the medium can be observed as 
a modification to the jet fragmentation function in heavy 
ion collisions. The fragmentation function is defined as 
D{z) = iv^^^, where z = p'^/j^"^; pj<=' is the initial 

jet momentum, and p'^ is the momentum of a hadronic 
jet fragment. Experimentally, this is accessible using di- 
rect photon-hadron correlations, where p^ « p!^^. This 
balance is only approximate due to the transverse mo- 
mentum, /ct, of the colliding partons inside nucleons, 
which on average introduces a transverse momentum im- 
balance and acoplanarity to the photon and its opposing 

jet [nlliilil. 

Several energy loss models [13, (HI predicting direct 
photon-hadron correlations only track the medium in- 
duced parton splitting of the leading parton. Other 
models follow the lost energy, leading to an increase in 
low momentum (soft) particle production. In particular, 
Borghini and Weidemann 16] use the modified leading 
log approximation (BW-MLLA) and local parton hadron 
duality to first reproduce the measured fragmentation 
function in e+e~ data. Modeling the energy loss in the 
medium as an increased parton splitting probability, they 
calculate the suppression of high jet fragments, as well 
as the redistribution of energy to lower px fragments and 
resulting enhancement at low z. The resulting Raa re- 
produces the PHENIX tt'^ measurement for 0-10% cen- 
tral events. The yet-another-jet-energy-loss model (Ya- 
JEM) traces the energy lost via gluon radiation and 
redistribution to soft particle production, predicting a 
suppression of particles at high z and an enhancement 
at low z. This calculation has been done specifically for 
"fdir-h, making it directly comparable to this data. The 
predicted Iow-zt enhancement has not yet been observed 
within the statistical and systematic limitations of pre- 



viously published data [18, 

In this letter, we report fragmentation functions mea- 
sured in Au-t-Au and p+p collisions determined from the 



yield of hadrons recoiling opposite to direct photons (i.e. 
the "away-side"). The extraction of a purely direct- 
photon sample is complicated by the presence of pho- 
tons from meson decays (dominantly tt^ — >■ 77), which 
must be removed from the inclusive photon-hadron corre- 
lations. PHENIX has previously established the extrac- 
tion of direct photon-hadron correlations via a statistical 
subtraction procedure in Au-f AuJlSf collisions and via 
an isolation cut in p+p collisions [11]. 

This analysis includes 3.9 billion minimum bias 
Au-I-Au events collected by PHENIX in 2007 and 2.9 
billion in 2010, after quality cuts. The p+p data set 
comprises 0.5 billion photon-triggered events collected 
in 2005 and 2006, corresponding to total recorded in- 
tegrated luminosities of 3.8 (2005) and 10.7 (2006) pb"\ 
respectively. Details on the p+p measurement were pre- 
viously presented in ll|. The kinematic reach and im- 
proved statistical precision of both data sets allow us to 
extend previous measurements 18, [l9| . reaching a low 
momentum fraction, z w 0.1, where interplay between 
the medium and the deposited energy may be impor- 
tant \7]. 

The Au-f Au minimum bias events are triggered by par- 
ticles firing the beam-beam counters (BBCs), which are 
arrays of Cerenkov counters covering 3.1 < \r]\ < 3.9 
and 2tt in azimuth. These BBCs are also used to de- 
termine the collision centrality and the collision vertex 
position along the beam direction. The 0-40% most cen- 
tral collisions are presented here. Photons and hadrons 
are measured in two central spectrometers spanning 7r/ 2 
in azimuth and ±0.35 units of pseudorapidity each |20|. 
The photons are measured in one of two electromagnetic 
calorimeters 2l| and charged hadrons are measured by 
reconstructing tracks in the Drift Chambers and Pad 
Chambers [22]. 

In Au-|-Au, a statistical subtraction determines the di- 
rect (i.e. nondecay) photon-hadron correlations from the 
measured inclusive photon-hadron correlations. Using 
the measured associated hadron yield per inclusive pho- 
ton, line = l/NincdN'^~^'-"'' /dA(j), and per decay photon, 
Ydoc, the associated yield per direct photon, Ydir, is de- 
termined by [l8| : 



Yaw 



R'y 1 



(1) 



Here Rj is the ratio of inclusive photons to decay pho- 
tons, reported by PHENIX in 0. 

Inclusive photon-hadron correlations are determined 
from the distribution of photon-hadron pairs as a func- 
tion of their azimuthal angular separation, Ac/). The dis- 
tribution of real pairs is divided by photon-hadron pairs 
in mixed events to correct for the PHENIX acceptance. 

The conditional, or per trigger, associated yield is ex- 
tracted after subtraction of photon-hadron pairs from the 
bulk underlying event 13|. Such particles are expected 



5 



to be correlated to one another only through the bulk 
anisotropy of the event, which is conventionally charac- 
terized by the Fourier coefficients f„, and are removed 



from the inclusive and _decay photon yields using the 
previously measured vi 



23[, and neglecting higher order 



terms, according to: 



J 



1 diVP^ 



1 r <;/^A0 _ ,^[i+2^„* a^,,3(2^^)]\ 



(2) 



r 



where the subscripts t and a refer to trigger and associ- 
ated particle, e° is the detection efficiency for the asso- 
ciated particle, and bo indicates the level of background 
pairs. The {V2V2) term modulates the background rate 
as a function of A4>, to account for correlations arising 
from flow of the bulk. 

The potential effect of ignoring in this extraction 
was studied by extrapolating the PHENIX hadron V3 
measurements (23l | . Including a modulation of the back- 
ground by this V3 in addition to the V2 results in a change 
in the away-side yield on the order of a few percent, de- 
pending on A(j) and px- The resulting background shape 
uncertainty is minimal for the highest hadron px selec- 
tions used here due to the low level of combinatorial 
background. The V3 effect is included as an additional 
systematic uncertainty on the background subtraction. 
An absolute background normalization is used to fix the 
background level, bo, as described in 24]. A GEANT sim- 
ulation of the detector determined the acceptance and 
efficiency for the measured charged hadrons, e". The 
uncertainty on e° leads to an 8.8% overall and 8.0% nor- 
malization uncertainty on the yields for the Au-|-Au and 
p+p data, respectively. 

To measure the decay photon contribution, n'^-h cor- 
relations are constructed following the same method as 
above using tt^s as the trigger. The tt^s are reconstructed 
from photon pairs whose invariant mass is within the 
window of 0.12-0.16 GeV/c^. The n'^-h correlations are 
translated into decay photon-hadron correlations accord- 
ing to a Monte Carlo study of the probability that a 7r° 
with a given px produces a decay photon within a certain 
Pt bin. This procedure is explained in detail in 18 1. 

In p+p collisions, jdir-h yields were measured using 
an isolation cut, and removing decay photons from the 
inclusive sample on an event- by- event basis ll|. In the 
analysis, photons which combine with another photon in 
the event to produce a mass within the tt" or rj mass 
windows were rejected. Next, an isolation cut was ap- 
plied, requiring that the transverse electromagnetic en- 
ergy and charged track momentum within a cone of 0.3 
rad around the photon be less than 10% of the photon 
energy. Finally, a statistical subtraction similar to that in 
the Au-t-Au analysis eliminated contributions from decay 
photons which appeared isolated or whose decay partner 



was lost due to finite detector acceptance or efficiency. 
In p+p collisions, the underlying event is subtracted as- 
suming the yield of photon-hadron coincidences is zero 
at the minimum point in the correlation as a function of 
A(f> (ZYAM). The lowest three points outside the isola- 
tion cut region are averaged; as there is no flow in p+p 
collisions, the background is assumed to be flat in A(f>. 
By first eliminating photons from other sources event-by- 
event, the signal to background ratio is improved and fi- 
nal uncertainties are reduced. The p+p results using the 
isolation cut agree with a statistical subtraction anal- 
ysis in p+p, but have smaller uncertainties. The high 
multiplicity of the underlying event makes it difficult to 
perform an isolation cut in Au-|-Au, so a statistical sub- 
traction procedure is used instead. 
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FIG. 1: (Color online) A(j) distribution for various associated 
^ bins for direct photons (circles) for the 0-40% most central 
Au-f Au collisions and the p+p reference (squares) in all pan- 
els. Panels (b), (d), and (f) are multiplied by a factor of 2 as 
indicated. 
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FIG. 2: (Color online) The top panel shows per trigger yield 
as a function of ^ for p+p collisions (squares) and 0-40% most 
central Au+Au collisions (circles). The points are shifted for 
clarity. For reference, the dependence on zt is also indicated. 
The bottom panel shows 7aa, the ratio of Au+Au to p+p 
fragmentation functions. Also shown are predictions from 
BW-MLLA tlQl (dashed line), calculated at Ejet = 7 GeV 
with fmed = 0.8 selected for 0-10% central Au+Au and from 
YaJEM-DE [H, ^ (dot-dashed curve) for 0-40% centrality 
and trigger photons from 9-12 GeV/c, both for the full away- 
side {\A(j) -n\ < n/2). 



In order to study the jet fragmentation function, D(z), 
associated hadron yields are determined as a function of 
Zt — Pt/Pt' ^^'^ ratio of the associated hadron trans- 
verse momentum, p^, to the trigger photon transverse 
momentum, p^. Here zt ~ z, since direct photon trig- 
gers balance the opposing jet. To focus on the low zt 
region, one can express the fragmentation function as a 
function of the variable, ^ = ln{l/zT)- To extend the 
accessible zt range, hadrons from 0.5 < px < 7.0 GeV/c 
are used in combination with a single 5 < p^ < 9 GeV/c 
photon bin. 

Figure [T] shows azimuthal pair angle distributions for 
the extracted direct 7-/1 correlations in 0-40% central 
Au+Au collisions as well as comparison with the direct 
7-/1 correlations in p+p. Unlike on the away-side, on the 
trigger side {\A<j>\ < 7r/2) the direct 7-/1 correlations in 
Au+Au show a negligible yield, indicating that the sta- 
tistical subtraction method indeed yields direct photons 
and that the yield of fragmentation photons in Au+Au 
is negligible within uncertainties. 



On the away side the associated particle yield is vis- 
ible, and there is significant variation when comparing 
the correlations in Au+Au to p+p. To further quan- 
tify this variation, the yields are integrated over A(f) for 
Itt — A(f)\ < 7r/2, as a function of ^, to obtain the effective 
fragmentation function. The top panel of Fig. [5] shows 
the integrated away-side yields in Au+Au and p+p as cir- 
cles and squares, respectively. The statistical error bars 
include the point-to-point uncorrelated systematic uncer- 
tainty from the background subtraction, while the boxes 
around the points show the correlated uncertainties. For 
reference, the dependence on zt is also indicated as the 
upper scale axis label. 

To study medium modification of the jet fragmenta- 
tion function, we take a ratio of the ^ distribution in 
Au+Au to p+p. This ratio, known as Iaa, is shown 
in the bottom panel of Fig. [2] and can be written as 
Iaa = yAu+Au^ yp+p^ y^^^^i of the global scale uncer- 
tainty cancels in this ratio, but there is a remaining 6% 
uncertainty. In the absence of modification, Iaa would 
equal 1. The data instead indicate suppression at low ^ 
and enhancement at higher ^. Including all systematic 
uncertainties the x^/dof value for the highest 4 points 
compared to the hypothesis that Iaa = 1 is 17.6/4, cor- 
responding to a probability that Iaa is 1.0 for ^ > 0.8 of 
less than 0.1%. 

The dashed curve in the bottom panel of Fig. [5] shows 
Iaa calculated at Ejet = 7 GeV using the BW-MLLA 
model in medium and in vacuum. The vacuum calcula- 
tion agrees well with the measured ^ distribution in e+e~ , 
and the in-medium conditions reproduce the measured 7r° 
Raa at high-pT for 0-10% central Au+Au events [l^ . 
The dot-dashed curve shows Iaa predicted by YaJEM- 
DE 25| for trigger photons from 9-12 GeV/c for the same 
centrality range (0-40%) as the present data [2^. Both 
models, which include all away-side jet fragments, show 
suppression at low ^ due to parton energy loss in Au+Au 
collisions, and increasing Iaa with increasing ^. In both 
cases, this is due to the lost energy being redistributed 
into enhanced production of lower momentum particles. 

The suppression of Iaa at high zt and enhancement 
at low Zt seen in these models agrees with the quali- 
tative trend in the data. However, the models do not 
reproduce the location in ^ where transition from sup- 
pression to enhancement is observed. Understanding the 
details of this transition can lead to better understand- 
ing of how lost energy is being redistributed. One such 
detail is how Iaa depends on the angular distribution of 
particles about the away-side jet axis. The top panel of 
Fig. [3] shows Iaa in three integration ranges. Reducing 
the integration range from |A0 — 7t\ < tt/2 reduces the 
observed enhancement and shifts the effect to higher ^. 
If the integration range is restricted to \A(j) — n\ < tt/6, 
the enhancement for ^ > 1.0 becomes negligible, while 
still showing significant suppression for ^ < 0.8. 

To better quantify the angular range of the enhance- 
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FIG. 3: (Color online) The top panel shows the Jaa for the 
full away-side (|A(j!) — 7r| < tt/2) (circles) and for two restricted 
away-side integration ranges, \A(j} — 7r| < n/3 (squares) and 
\A(j) — 7r| < 7r/6 (triangles). The points are shifted for clarity. 
The bottom panel shows the ratio of the Jaa for \A<j} — tt] < 
n/2 to \ A4i - tt] < 7r/6. 



ment, we can look at the ratio of /aa's with different 
integration ranges, where some of the statistical and sys- 
tematic uncertainties common to all Iaa cancel. The 
bottom panel of Fig. [3] shows the ratio of the full away- 
side integration range to the \A(j> — 7r| < tt/6 case. From 
this ratio it is clear that there is a significant variation 
in observed Iaa as a function of the integration range. 
The average ratio for ^ > 0.8 is 1.9±0.3(stat) ±0.3(syst), 
indicating that the enhancement in Iaa seen at large ^ 
is predominately at large angles (|A^ — tt\ > 7r/6). 

In summary, we have presented evidence for medium 
modification of jet fragmentation, measured via compar- 
ison of direct photon-hadron correlations in ^s„„ = 200 
GeV Au-l-Au and p+p collisions. The ratio of Au-f-Au to 
p+p yields indicates that particles are depleted at low ^ or 
high momentum fraction, zt, due to energy loss of quarks 
traversing the medium. The ratio exhibits an increasing 
trend toward high ^, exceeding one at ^ > 1.0. Restrict- 
ing the away-side azimuthal integration range reduces the 
enhancement at high ^ significantly. This suggests that 
the medium enhances production of soft particles in par- 
ton fragmentation, relative to p+p, preferentially at large 
angles. 
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